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Six Poisson autoregressive models of order p [PAR(p)] of daily wildland arson ignition counts are 
estimated for five locations in Florida (1994-2001). In addition, a fixed effects time-series Poisson model 
of annual arson counts is estimated for all Florida counties (1995-2001). PAR(p) model estimates reveal 
highly significant arson ignition autocorrelation, lasting up to eleven days, in addition to seasonality 
and links to law enforcement, wildland management, historical fire, and weather. The annual fixed 
effects model replicates many findings of the daily models but also detects the influence of wages and 
poverty on arson, in ways expected from theory. All findings support an economic model of crime. 
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Arsonists set i^Tnillion fires each year in the 
United States, resulting in over $3 billion in 
damages, 500 fatalities, and thousands of in¬ 
juries (TriData Corporation). From the stand¬ 
point of wildfire, arson comprises a significant 
share of all wildland fires in many parts of 
the country, especially in populated regions. 
Hall reports that about 10% of property lost 
to fire is attributable to outdoor firesetting. 
Arson wildfires adversely affect wildland man¬ 
agement of timber, water, recreation, grazing, 
and biodiversity. 

Wildland arson threatens lives and the 
stability of the communities that depend on 
land-based social goods and services. The 
arson-ignited 2002 Hayman fire, near Denver, 
burned 138,000 acres, destroyed over 100 res¬ 
idences, and generated damages and costs for 
governments that totaled over $100 million. 
Evacuations and other disruptions created fur¬ 
ther economic losses and transfers that totaled 
into the tens of millions of dollars (Kent et al.). 
Although research is limited on the causal fac¬ 
tors for wildland arson, preliminary evidence 
suggests that arson wildfires may cause a dis- 
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proportionate share of losses and costs because 
they are more commonly ignited near values at 
risk (Butry, Pye, and Prestemon). 

In spite of its potential importance, research 
into the factors influencing wildland arson has 
been limited to only a few published stud¬ 
ies (e.g., Donoghue and Main). This contrasts 
with a more abundant and recent literature 
on property crime. For example, rural and ur¬ 
ban property crime has been linked to eco¬ 
nomic conditions, including poverty (Arthur; 
Brotman and Fox; Hannon; Hershbarger and 
Miller; Neustrom and Norton), unemploy¬ 
ment and wages (Carlson and Michalowski; 
Burdett, Lagos, and Wright; Gould, Weinberg, 
and Mustard; Spillman and Zak), and race 
(Viscusi). A new body of research has 
connected crime to law enforcement effort 
(Cameron; Corman and Mocan; Di Telia and 
Schargrodsky; Eck and Maguire; Fisher and 
Nagin; Marvell and Moody). 

Potentially more important than the limited 
research linking arson to hypothesized causes, 
wildland arson has not been analyzed as a 
temporal process. This is in spite of colloquial 
knowledge of its temporal clustering and in 
spite of research that has identified cluster¬ 
ing and seasonality for other kinds of crimes 
(e.g., Farrell and Pease; Johnson and Bowers; 
Bowers and Johnson; Surrette; Townsley and 
Pease; Townsley, Homel, and Chaseling). One 
reason for a lack of progress in capturing 
temporal clustering is that only recently have 
valid statistical approaches been developed 
that can account for it in event models (e.g., 
Chang, Kavvas, and Delleur; Zeger; Harvey 
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and Fernandes; Cameron and Trivedi; Brandt 
et al.; Brandt and Williams 2001). 

The correct specification of an event count 
model (see Cameron and Trivedi) that includes 
an underlying autoregressive data generation 
process for daily wildland arson ignitions 
would be justified for statistical, as well as be¬ 
havioral, reasons (see Brandt and Williams 
2001). Statistical consistency, unbiasedness, 
and efficiency objectives of count modeling 
require explicit accounting for autoregressiv- 
ity (Brandt and Williams 2001). Structure ar¬ 
sonists often set multiple fires in a short time 
frame (Sapp et al.), leading to a hypothe¬ 
sis that wildland arsonists demonstrate the 
same tendency. Also, youth-caused property 
crime has been shown to include “copycat” el¬ 
ements (Surrette), so wildland arsonists might 
be prone to this behavior, as well. Augment¬ 
ing both serial and copycat phenomena is the 
informational content of recent arson fires: ar¬ 
sonists could view a successful wildland ar¬ 
son ignition as a signal that environmental 
conditions are favorable for successfully start¬ 
ing new fires, a plausible encouragement for 
more. 1 

The primary goal of this research is to 
model wildland arson to account for tempo¬ 
ral clustering and in order to test theories in 
both criminology and wildland management. 
We estimate two kinds of count data mod¬ 
els using county-level data in Florida. The 
first is a Poisson autoregressive model of or¬ 
der p (Brandt et al.; Brandt and Williams 
2001) of daily wildland arson ignitions. This 
model tests the hypothesis that arson is tem¬ 
porally clustered into multiple-day outbreaks. 
The second is a fixed-effects panel Poisson 
model of annual wildland arson ignitions for 
all Florida counties. Although this model can¬ 
not capture fine time-scale clustering, it intro¬ 
duces annual variation in wildland arson that 
may better account for how arson is affected 
by variables that change slowly or that can¬ 
not be expressed accurately at the daily time 
scale. 


1 Alternatively, autoregressive patterns in wildland arson are 
only due to environmental factors that are not observable by 
the analyst or that capture imprecision at finer spatial or tem¬ 
poral scales. In other words, arson ignition attempts could occur 
at a constant rate throughout a year, but it is only during cer¬ 
tain weather that ignitions are successful. Lagged ignitions might 
then proxy for the unobservable local weather. If this is the case, 
then our modeling addresses successful arson ignitions, not arson 
attempts. 
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Conceptual Framework 

Becker specified person fs decision on the 
commission of a crime as 

(1) o, = 0/K, //, m) 

where tt,- is the probability of being caught and 
convicted, /,• is an income-equivalent loss ex¬ 
perienced by the offender for being caught and 
convicted, and measures all the other in¬ 
fluences on the offense. Equation (1) can be 
called a crime function (Fisher and Nagin). The 
first derivatives of O, with respect to tt, and /, 
are negative. 

Now, define the arsonist’s psychic and in¬ 
come benefits from illegal firesetting as g, and 
the production cost for the firesetting as c,. 
As described by Burdett, Lagos, and Wright, 
and ignoring issues of risk aversion, the loss 
from being caught and convicted of commit¬ 
ting the crime, /,, is a positive function of in¬ 
come while employed, so that/, =/, (w,, W,), 
where w, is the wage rate (Gould, Weinberg, 
and Mustard) and W, is the employment status. 
Adapting Becker (footnote 16), the prospec¬ 
tive arsonist's expected utility from success¬ 
fully starting a wildland arson fire may be 
expressed as 

(2) EUi(Oj) = tt jUi(gi - Ci - /(IT, w)) 

+ (1 - TT,)f/,(g,- - a). 

As wages rise, for example, the expected net 
utility from arson declines, lowering the prob¬ 
ability that an arson fire will be set: 3 EUi (O,)/ 
dwi = , n(dEUi/df j)(df j/dwi) < 0. 2 

As with opportunity costs, the production 
costs of firesetting can be more elaborately 
described. For example, c, may be a func¬ 
tion of time available (Jacob and Lefgren); 
fuels and weather (Gill et al., Vega Gar¬ 
cia et al., Prestemon et al.); whether the 
person is unemployed (creating a lower op¬ 
portunity cost of crime); and information 


2 The marginal utility of committing arson would be larger if the 
opportunity cost of work were included as an additional cost. As¬ 
suming that Ci = Ci(wi), and 3 a (wi)/dwi > 0, then dEUi ( Oi)/dwi 
= tt (dEUi/dfi)(dfi/dwi) + ( dEUi/dCi)(dCi/dwi ) < 0. Because the 
second term on the right-hand side of this expression is also nega¬ 
tive, the marginal effect of a wage increase on utility is even larger 
than if the opportunity cost of time to commit the arson crime 
is 0. 
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on other successfully ignited arson wildfires, 
which lowers the production cost by rais¬ 
ing the success rate. 3 Anything that raises 
the crime production cost, lowers the ex¬ 
pected utility of the crime: 3 EU-, (() t )/ 
da = Tr(dUi/dCi ) + (1 — tt )(3[/,'/3 c ,) < 0. 

Similarly, following Burdett, Lagos, and 
Wright, tt can be expressed as a function of 
law enforcement effort. As Fisher and Nagin 
pointed out, aggregate crime may be simul¬ 
taneously determined with law enforcement, 
so not accounting for this simultaneity when 
it is present can cause a positive bias in a 
measured statistical effect of law enforcement 
on crime (Cameron; Marvell and Moody; 
Eck and Maguire). Becker describes how 
law enforcement effort (spending) to reduce 
crime is jointly determined with the crime 
rate. Recent research has used instrumental 
variables methods to resolve this problem in 
estimating crime functions. In any case, it is 
likely that bureaucratic delays would mean 
that law enforcement agencies would find it 
difficult to quickly change enforcement rates 
in response to higher crime rates. There is a lag 
between hiring new officers and effectiveness 
of the new officers in the field (Corman and 
Mocan). 4 Gould, Weinberg, and Mustard 
showed that the kind of endogeneity involved 
with one measure of law enforcement (police 
presence or spending) is small. We contend 
that aggregate law enforcement is unlikely to 
change in response to wildland arson, given 
that most crime is not wildland arson. We 
therefore believe that assuming exogeneity is 
justified. 5 


3 Alternatively, Burdett, Lagos, and Wright describe -tt as the 
probability of being caught, conditional on having an opportunity 
to commit a crime. Following this, weather and fuels could provide 
the opportunity, so that tt should be a function of law enforcement, 
weather, and fuels. In our empirical specification, this difference is 
not important. 

4 An ideal model would instrument E before including it in the 
model as a contemporaneous predicted variable, but the avail¬ 
able “best” instruments for E at the county level or at the two- 
county level using daily data are limited or unavailable. We also at¬ 
tempted to control for simultaneity by directly including an instru¬ 
ment in our annual model—a measure of so-called “index crime” 
(Florida Department of Law Enforcement 2004)—and found no 
large effects on parameter estimates or significances. Although 
Gould, Weinberg, and Mustard did not find that the bias affected 
significance or direction of effect of many other predictors of 
crime, some doubt remains and suggests an avenue for additional 
research. 

5 In pooled time-series cross-sectional data on annual arson ig¬ 
nitions for all 67 Florida counties, 1995-2001, wildland arson ig¬ 
nitions per capita and police per capita are negatively correlated 
(r = —0.14). This contrasts with the positive correlation between 
the crime index (Florida Department of Law Enforcement 2004) 
and police per capita (+0.29), providing some evidence of joint 


Empirical Specifications 

A PAR(p) Model of Daily Wildland 
Arson Ignitions 

The PAR(p) model, generalizing work by 
Grunwald, Hamza, and Hyndman, was origi¬ 
nally developed by Brandt and Williams (2001) 
to model persistence in annual patterns of 
presidential vetoes and purse snatchings. In 
our application to daily data, we start by sum¬ 
ming the daily decisions on crime made by all 
persons (i — 1 to /) in location j. The aggre¬ 
gate arson fire outcome of these I decisions on 
day t in location j is a count of arson ignitions, 
yjj. The PAR(p) model hypothesizes that the 
observed count is drawn from a Poisson distri¬ 
bution conditional on nij t , 

e~ m t' 1 

(4) Pr[y;,, | »»;,,]= JJ , 

yj.t- 

where nij t — E[y ( , | Yy t \ is the conditional 
mean of a linear AR(p) process. The expected 
count can be described as 


E hit I Y i,t~ i] = X! p DyD-i 

i =1 

1 - E P/,z)exp(x , ; ,(i / ) 


where Xjj is a vector of independent variables 
(including a constant), (J, is a vector of associ¬ 
ated parameters, and the p / ( ’s are the autore¬ 
gressive parameters. 

The density of (5) has a gamma-distributed 
conjugate prior such that 

(6) Pr(my, t 110,,_i) = T(oj;,.pn;,,-!, 

where my r _i and the variance cr?, are both 
positive, F(-) is the gamma distribution, and 
myt -1 = E[y/,r I Yj,t-\\ and CT? t _ x = V\yy t \ 
The likelihood equation associated 


determination, and the negative correlation between the crime in¬ 
dex and wildland arson per capita (—0.19), which is circumstan¬ 
tial evidence that law enforcement is exogenous to wildland arson 
ignitions. 
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with this model is (omitting for simplicity the 
location subscript /): 

(7) 

°r-i | Y t - 1) 

T T 

= In Y\ Pr (>v I Yt) = X! lnr ( fJ ^r-! m r-1 + Jr) 

r=l f=l 

- r (yr +1) - r(o-, 2 _ 1 m,_i) 

+ cr}_ 1 m t -1 ln((rp_^) 

+ y ? )ln(l + (T, 2 .!). 

The associated instantaneous (short-run) im¬ 
pact (marginal effect) of a change in variable 
Xk in x on the count is described as 


(8a) 


dm j it 
dx j,k,t 


p 

£p, 


i =1 


cxpfx'jijp 


'j,k- 


The elasticity of the ignition count, at some 
reference point in the relationship between 
the count, and the variable in question corre¬ 
sponding with the impact, could be calculated 
simply by multiplying (8a) by Here, 

Xj i; and ni j could be the mean values for Xjj ( 
and trij observed in the data. The associated 
long-run impact of a change in variable in x 
on the count in location j is described as 


(8b) 


dm j,t / dx 

( 1 -Ef„p w ) 




and the elasticity would be calculable in a man¬ 
ner analogous to the instantaneous effect, i.e., 
cxp(xj ( P ; ) $j,kXj,k /nij • Note, however, that the 
expected values of the estimated parameters in 
the PAR(p) model do not necessarily equal the 
expected values of the estimated correspond¬ 
ing parameters in the Poisson (unless all p / ( - = 
0). Brandt and Williams (2001) illustrate, per a 
Monte Carlo simulation, that if the process is 
truly a PAR(p) process, then a Poisson model 
will underestimate the true long-run impact. 

Our daily arson ignition model is esti¬ 
mated for five high arson locations in Florida. 
These five are comprised of four two-county 
pairs (Duval-St. Johns, Flagler-Volusia, Taylor- 
Dixie, Sarasota-Charlotte) and one single 
county (Santa Rosa). The State of Florida iden¬ 
tifies these as five homogeneous areas having 
the largest concentration of arson in the state. 
As long as our assumptions that behavioral or 
underlying relationships between dependent 
and independent variables hold statewide, this 


sample would be representative of aggregate 
arson behavior in Florida. Otherwise, our find¬ 
ings are strictly applicable just to the five loca¬ 
tions. 6 

The full specification of a daily arson crime 
function, consistent with Becker; Burdett, 
Lagos, and Wright; Gould, Weinberg, and 
Mustard; and Jacob and Lefgren, includes 
variables measuring law enforcement effort, 
economic conditions, fuels and fuels manage¬ 
ment, weather, and daily and seasonal factors. 
Hence, our model includes the location’s daily 
interpolation of law enforcement effort, E 
(i.e., tt, = it (£))); a daily interpolation of the 
percent of households below the poverty line 
(i.e., the poverty rate), an index of inequality 
or social disorganization, which also captures 
the relative cost of crime; a daily interpolation 
of the location’s unemployment rate (aver¬ 
aged over the months of the year); a daily 
interpolation of the state’s average annual 
retail wage rate; the running total extent 
(acres) of lagged wildfire in the previous zero 
to two and three to six years; the running total 
permitted prescribed burning extent (acres) in 
the previous zero to one and one to two years; 
a daily observation of a drought measure, 
the Keetch-Byram Drought Index (KBDI); 
dummy variables that index weekend days and 
holidays (representing more time available for 
firesetting, possibly measuring lower opportu¬ 
nity costs of time); and dummy variables that 
index months, accounting for other sources of 
intra-annual seasonality. Further, the specifi¬ 
cation includes a time index that controls for an 
unknown set of slowly changing factors, such 
as population, the amount of wildland, and law 
enforcement technology. Finally, autoregres- 
sivity in wildland arson is captured by lagged 
arson wildfire ignitions. 7 The order of autore¬ 
gression varies by location, depending on the 
availability of sufficiently long-lasting arson 
clusters to permit parameter identification. 8 


6 The annual arson ignition model is how we seek to correct for 
this inferential limitation. 

7 We included total population in an initial empirical specifica¬ 
tion, but it was omitted due to nonconvergence. Hence, the primary 
effect of population (its absolute level) is contained in the inter¬ 
cept of our daily models. Population is included, however, in the 
annual model. 

8 GAUSS 5.0 with the Maximum Likelihood 5.0 Application 
Module (Aptech Systems, Inc.) was used to estimate the PAR (p), 
model. Parts of the programming code (designed for GAUSS 3.0) 
were provided by Brandt and Williams (2002), although we up¬ 
dated this code to make it compatible with the later version of 
GAUSS and our particular modeling framework. NLOGIT 3.0 
(Greene) was used to estimate the fixed effects Poisson model. 
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A pooled version of the PAR(p) individual- 
location model is also estimated for the five 
locations. The pooled version is estimated to 
test inferences with a larger dataset and to 
quantify more precisely the average amount of 
autoregressivity present. One way to estimate 
a pooled version would be to incorporate 
cross-sectional dummy variables, a PAR(p) 
equivalent to a dummy variable or a fixed 
effects panel time-series model. However, 
no statistical methods for the panel PAR(p) 
model are available. Here, we assume that 
parameters of individual-location versions of 
equation (5) differ across cross-sections only 
by a factor proportional to the population of 
the location. Multiplying the nonautoregres- 
sive terms on the right-hand side of (5) by 
population effectively controls for this pro¬ 
portional difference. Because of these interac¬ 
tions, the right-hand side of (5) also includes 
population on day t. The intercept is still re¬ 
quired for statistical consistency. The time pe¬ 
riod of inference varies by location but roughly 
covers January 8,1994 to December 31, 2001. 

Annual Arson Ignitions 

The annual model of arson ignitions is a fixed- 
effects panel Poisson process. In each year, the 
arson ignition count is the sum of the arson 
fire outcomes resulting from I daily decisions 
on firesetting across all days of the year. For 
cross-sectional units j = 1 ,...,/ and years t = 
1 ,,T, y,.T■ the density of this annual count 
is described as (Greene) 


where 


(12) —exp(oq) 


E 


V^T 

z2t=i yj, t 
t=i exp(x; T p) 


As in the case of the daily model of wildland 
arson, the annual specification includes vari¬ 
ables related to law enforcement, economic 
conditions, fuels and fuels management, and 
weather. Variables included are: current year 
police officers per capita, the poverty rate, the 
average annual unemployment rate, the state- 
level average annual retail wage rate, twelve 
separate years of lags of wildfire area burned 
(acres), the current and previous year’s per¬ 
mitted hazard reduction prescribed burning 
(acres), and the total pulpwood harvest vol¬ 
ume in each of the three previous years in 
the county (thousands of cubic feet). Average 
weather statewide is quantified by the year’s 
average Nino-3 sea surface temperature (SST) 
anomaly (°C), which is a measure of the El 
Nino-Southern Oscillation, and a dummy vari¬ 
able that accounts for the unusual El Nino 
of 1997-8 (see Prestemon et al.). To control 
for other trends related to population and ag¬ 
gregate urban growth (which also consumes 
fuels), we include county population and a 
state-level time trend. Finally, county fixed ef¬ 
fects measure cross-sectionally varying factors 
that did not change significantly over the pe¬ 
riod of inference, 1995-2001, such as ecological 
variables. 


Data 


( 9 ) f(yj, x \ Xj,l,Xj, 2 , ...,x yVr ) 

= i'(X/.TA' / . T P + ot ; ) 

where the ay’s are cross-sectionally related 
fixed intercepts. The fixed-effects Poisson 
count model 9 has an expected count, 

(10) E(>., t IX T ) = \ ; , T = exp(x' ; T p + ay). 
The likelihood equation to maximize is 

/ T 

(11) logL = EE —exp(ay)exp(x' T p) 

7=1 T=1 

+ ( x y, T P + «;) — lo g>’7,T ! 


Wildfire and prescribed fire permit data were 
obtained directly from the Florida Division of 
Forestry. Arson wildfires were those deemed 
by the Division as likely arson, but uncer¬ 
tainty means that an unknown number of 
fires were misclassified. 10 Local population es¬ 
timates were from the Florida Bureau of Eco¬ 
nomic and Business Research, while annual 
poverty data were from the United States 
Department of Commerce, Census Bureau. 
The Florida Department of Law Enforcement 
(2002) provided data on the mid-year count 
of full-time equivalent police officers in each 
county. The retail wage rate in our models was 
the state-level average for the year, from the 
United States Department of Labor (2004). 


9 We confine our discussion to the Poisson model. In empiri- 
cal estimates for this analysis, the standard alternative model, the 
negative binomial, which accounts for overdispersion, did not con¬ 
verge in estimation. 


10 Any wildfire of suspicious origin is investigated, and correct 
classification of its cause is highly likely, according to Division per¬ 
sonnel. Nevertheless, classification errors would result in an un¬ 
known degree of statistical inconsistency in parameter estimates. 
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Table 1. Summary Statistics for Five (Two- and Single-County) Locations in Florida, for the 
Daily Model 


Variables 



County Aggregate 



Dixie + 

Taylor 

Duval + 

St. Johns 

Sarasota + 
Charlotte 

Santa 

Rosa 

Volusia + 
Flagler 

Arson ignitions/day 






Mean 

0.18 

0.20 

0.25 

0.24 

0.20 

Maximum 

18 

8 

7 

17 

5 

Minimum 

0 

0 

0 

0 

0 

SD 

0.89 

0.58 

0.77 

0.88 

0.57 

Police officers 






Mean 

48 

1,738 

405 

654 

1,113 

Maximum 

50 

1,859 

445 

730 

1,217 

Minimum 

45 

1,564 

355 

568 

1,005 

SD 

1 

84 

22 

46 

60 

Poverty rate (% of population) 





Mean 

20.09 

15.04 

8.82 

11.25 

13.78 

Maximum 

23.06 

16.73 

9.78 

12.50 

14.87 

Minimum 

15.52 

14.02 

7.83 

4.89 

12.02 

SD 

2.50 

0.82 

0.62 

1.37 

1.02 

Unemployment rate (% of population) 





Mean 

7.94 

3.50 

2.87 

3.85 

3.68 

Maximum 

11.75 

4.97 

4.33 

5.70 

5.24 

Minimum 

5.26 

2.53 

1.85 

2.81 

2.64 

SD 

1.40 

0.46 

0.54 

0.50 

0.62 

Real wage rate (2001 $/yr) 





Mean 

16,564 

16,745 

16,719 

16,765 

16,723 

Maximum 

17,468 

17,669 

17,631 

17,707 

17,631 

Minimum 

16,146 

16,146 

16,146 

16,146 

16,146 

SD 

375 

480 

464 

496 

464 

Wildfire (ac/yr) 






Mean 

4,126 

4,631 

5,478 

2,663 

42,159 

Maximum 

21,605 

22,661 

12,836 

5,082 

263,026 

Minimum 

70 

411 

843 

644 

754 

SD 

7,559 

6,836 

3,266 

1,438 

91,688 

Hazardous burn permits (ac/yr) 





Mean 

3,056 

1,348 

5,034 

14,236 

4,203 

Maximum 

9,055 

4,676 

20,316 

38,061 

14,719 

Minimum 

24 

0 

0 

571 

51 

SD 

2,366 

1,272 

6,125 

8,986 

3,968 

Population (thousands) 





Mean 

18.3 

785.7 

338.4 

109.9 

436.0 

Maximum 

18.9 

822.3 

356.7 

122.7 

461.6 

Minimum 

17.6 

745.3 

322.3 

96.5 

412.2 

SD 

0.4 

23.2 

10.0 

8.4 

15.1 

Keetch-Byram drought index 





Mean 

314 

295 

431 

222 

312 

Maximum 

721 

733 

783 

681 

694 

Minimum 

0 

1 

4 

0 

1 

SD 

217 

166 

195 

183 

181 


County unemployment data were from the 
United States Department of Labor (2002). 
The current day’s KBDI was constructed using 
an algorithm (Keetch and Byram) from rep¬ 
resentative weather station data in the study 
area, which were collected by the National 
Climatic Data Center and provided by Earth- 
Info, Inc. 


Table 1 presents summary statistics for 
nonseasonal, nontrend variables that we 
use in the PAR(p) modeling. Average daily 
arson ignitions are similar across the five 
locations during the sample period, although 
the maxima observed in each varied more—up 
to 18 in one case in the Dixie-Taylor county 
aggregate. The full-time equivalent police 
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officers per county vary from under 50 to 
over 1,700. Poverty rates in the Dixie-Taylor 
aggregate are double those observed in the 
Sarasota-Charlotte aggregate, and those two 
locations represent the highest and lowest 
rates of unemployment. The state-level wage 
rate does not vary across locations, of course, 
so differences in summary statistics for each 
location are caused by differences in sample 
periods, which are constrained by the avail¬ 
ability of weather data. Wildfire area burned 
varies by nearly two orders magnitude across 
locations (large fires in Volusia and Flagler 
counties in 1998 distort that location’s average 
and standard deviation upward), but pre¬ 
scribed fire rates are less variable. Population 
varies widely—Dixie and Taylor are mainly 
rural counties, while most of the rest are more 
urbanized, with higher pressures on the wild¬ 
land resource in the latter group. As is clear 
from the data on the KBDI, weather-related 
fire danger is highly variable over the sample 
periods for all locations, with values ranging 
from the index minimum of 0 (soil saturation) 
to near its maximum of 800 (maximum soil 
moisture deficiency or maximum drought 
conditions) (Keetch and Byram). 

Additional data not shown in table 1 but 
used in the annual arson fixed-effects model¬ 
ing include pulpwood harvest data, obtained 
by special request from the Southern Research 
Station of the United States Department of 
Agriculture, Forest Service, and the Nino-3 
SST anomaly, obtained from the National 
Oceanic and Atmospheric Administration. 
Summary statistics by county aggregate on 
these and the other variables included in 
the county-level annual fixed-effects model 
are available from the authors upon request. 
We note here that the data summarized in 
table 1 for the daily models approximate 
the statewide variability of most modeled 
variables in the annual models. 

Results 

Daily PAR(p) Model Results 

Daily time-series arson ignition models are 
broadly significant, and the included variables 
explain a large share of the variation in ig¬ 
nitions in each location (table 2) and in the 
pooled data model (table 3). Although PAR(p) 
model estimates were attempted with up to 
twelve autoregressive terms for all locations, 
parameter identification restrictions allowed 
only smaller versions to be estimated. The 


most commonly statistically significant deter¬ 
minants of arson ignitions are police officers 
per capita (negatively related), a dummy vari¬ 
able for Saturday (with the expected positive 
effect), the high-fire season months of the year 
(positively), and the autoregressive terms. The 
poverty rate, expected to have a positive influ¬ 
ence, is statistically significant in two cases—in 
the Volusia-Flagler estimate and in the pooled 
estimate, in both cases with the expected sign. 
Unemployment and retail wages are related to 
arson in directions not found in other research; 
they are usually statistically insignificant, but 
wages are unexpectedly and significantly pos¬ 
itively related to arson in the pooled model. 
The number of statistically significant autore¬ 
gressive terms varies between six and eleven, 
the most in the pooled model. Tests of model 
restrictions that the autoregressive parameters 
are jointly zero are rejected at less than the 1 % 
significance level in all daily models. 

Significance tests show that arson wildfire 
ignitions also vary with weather and fuels. 
Fligher Keetch-Byram drought indices are 
linked to higher arson counts, other variables 
held constant. The lagged wildfire and pre¬ 
scribed fire variables, both being inversely 
related to fuel loads, are only occasionally 
statistically significant. In those few cases, 
the measured effects of these variables are 
negative, as expected. 

To assess the net effect of explicitly ac¬ 
counting for autocorrelation in daily arson 
ignitions, we compare average elasticities 
produced by the individual-location and 
pooled PAR(p) models with counterpart 
nonautoregressive Poisson alternative speci¬ 
fications (table 4). When the parameter esti¬ 
mates of the PAR(p) model are statistically 
significant, usually so too are their nonautore¬ 
gressive counterparts. However, significance 
levels differ between the PAR(/;) and the Pois¬ 
son in several instances. For example, in the 
pooled specification, the Poisson model indi¬ 
cates that police are positively related to ar¬ 
son rates, while the PAR(p) model shows no 
significant influence (and a negative effect in 
two individual-location models). The Poisson 
model estimate provides evidence that wild¬ 
land arson ignitions are significantly and nega¬ 
tively linked to lagged wildfire and prescribed 
burning, while the PAR(p) model indicates no 
effect. These conflicting results suggest a re¬ 
lationship between some observed variables 
and the underlying autoregressive pattern, 
which the PAR(p) model specifically handles 
but the static model does not. In general, the 



Prestemon and Butry 


Wildland Arson as an Autoregressive Crime 763 


Table 2. Poisson Autoregressive Models of Maximum Order Estimable, Five Study Areas in 
Florida, Daily Counts of Wildland Arson Ignitions, 1994-2001 


Variables 


Model Locations 


Dixie- 

Taylor 

Duval-St. 

Johns 

Santa 

Rosa 

Sarasota- 

Charlotte 

Volusia- 

Flagler 

Constant 

52.69* 

-28.82 

-27.92** 

11.46 

-59.32 


(29.95) 

(85.24) 

(13.06) 

(30.24) 

(43.04) 

Police per capita 

-1.30** 

0.39 

0.15 

-4.16*** 

0.17 


(0.58) 

(2.82) 

(0.18) 

(1.20) 

(1.18) 

Poverty rate 

-0.16 

-0.02 

0.14 

2.96* 

0.95 


(0.62) 

(0.55) 

(0.20) 

(1.65) 

(0.77) 

Unemployment rate 

4.55 

9.34 

-3.57 

-7.15*** 

3.60 


(3.67) 

(11.24) 

(2.48) 

(1.70) 

(5.74) 

Real retail wage 

-1.65 

0.96 

1.12 

0.39 

2.21* 


(1.22) 

(1.19) 

(0.72) 

(0.74) 

(1-19) 

Wildfire years t to t — 2 

-0.45 

0.08 

-3.52*** 

-0.87 

-0.01 


(0.41) 

(0.63) 

(1.27) 

(0.53) 

(0.03) 

Wildfire years t — 3 to t — 5 

-0.75 

-0.50 

-1.29 

0.18 

-0.01 


(1.51) 

(0.59) 

(1.44) 

(0.24) 

(0.05) 

Prescribed fire year t 

3.49** 

2.07 

0.31 

-0.11 

1.54* 


(1.61) 

(2.18) 

(0.20) 

(0.18) 

(0.94) 

Prescribed fire year t — 1 

-0.25 

1.24 

-0.07 

-0.22 

0.04 


(1.01) 

(2.14) 

(0.16) 

(0.23) 

(0.66) 

KBDI 

0.39*** 

0.32*** 

0.26*** 

0.32*** 

0.30*** 


(0.10) 

(0.09) 

(0.08) 

(0.06) 

(0.09) 

Saturday 

1.18*** 

0.61** 

0.30 

0.060 

0.70** 


(0.28) 

(0.29) 

(0.28) 

(0.078) 

(0.31) 

Sunday 

0.46* 

0.11 

-0.02 

0.047 

0.27 


(0.27) 

(0.38) 

(0.25) 

(0.081) 

(0.32) 

Holiday 

0.36 

1.90*** 

-0.67 

0.10 

0.31 


(0.54) 

(0.31) 

(0.57) 

(0.16) 

(0.48) 

January 

0.33 

0.48 

2.06*** 

0.50 

0.47 


(0.48) 

(0.52) 

(0.52) 

(0.34) 

(0.49) 

February 

1.72*** 

1.40*** 

2.35*** 

0.46 

1.28*** 


(0.52) 

(0.49) 

(0.49) 

(0.34) 

(0.47) 

March 

1.50** 

1.17** 

1.40** 

0.61* 

0.44 


(0.61) 

(0.53) 

(0.55) 

(0.34) 

(0.56) 

April 

0.33 

1.51*** 

1.31** 

0.81** 

0.78 


(0.64) 

(0.50) 

(0.53) 

(0.34) 

(0.52) 

May 

0.24 

0.43 

0.79 

0.76** 

0.27 


(0.60) 

(0.55) 

(0.56) 

(0.36) 

(0.55) 

June 

0.50 

-0.43 

0.03 

-0.11 

-0.60 


(0.78) 

(0.81) 

(0.65) 

(0.54) 

(0.69) 

July 

-1.04 

-0.10 

-0.45 

-0.73 

0.35 


(0.82) 

(0.60) 

(0.61) 

(0.57) 

(0.61) 

August 

-0.33 

0.06 

-0.70 

-1.11* 

-0.75 


(0.87) 

(0.60) 

(0.62) 

(0.67) 

(0.68) 

September 

-0.90 

-1.26 

-0.46 

0.21 

-0.38 


(0.94) 

(1.07) 

(0.61) 

(0.48) 

(0.64) 

October 

-0.18 

0.48 

0.35 

0.36 

-1.28 


(0.75) 

(0.55) 

(0.57) 

(0.47) 

(0.82) 

November 

-0.87 

-1.00 

0.32 

-0.40 

0.35 


(0.72) 

(0.90) 

(0.57) 

(0.49) 

(0.55) 

Time 

0.25 

-0.12 

0.04 

0.39** 

0.05 


(0.28) 

(0.30) 

(0.05) 

(0.18) 

(0.18) 

Arson ignitions day t — 1 

0.009 

0.20*** 

q 29*** 

0.14*** 

0.22*** 


(0.038) 

(0.04) 

(0.05) 

(0.04) 

(0.05) 

Arson ignitions day t — 2 

0.036 

0.16*** 

0.038 

0.10*** 

0.091** 


(0.035) 

(0.04) 

(0.035) 

(0.04) 

(0.043) 
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Table 2. Continued 





Model Locations 




Dixie- 

Duval-St. 

Santa 

Sarasota- 

Volusia- 

Variables 

Taylor 

Johns 

Rosa 

Charlotte 

Flagler 

Arson ignitions day t — 3 

0.17*** 

0.055* 

0.060* 

0.055* 

0.067* 


(0.05) 

(0.033) 

(0.034) 

(0.032) 

(0.040) 

Arson ignitions day t — 4 

0.088** 

0.042 

0.13*** 

0.060* 

0.12*** 


(0.041) 

(0.030) 

(0.04) 

(0.033) 

(0.04) 

Arson ignitions day t — 5 

0.027 

0.013 

0.00 

0.051 

0.024 


(0.045) 

(0.029) 

(0.03) 

(0.034) 

(0.036) 

Arson ignitions day t — 6 

0.076* 

0.14*** 

0.043 

0.017 

on*** 


(0.046) 

(0.04) 

(0.030) 

(0.034) 

(0.04) 

Arson ignitions day t — 7 


0.087*** 

0.032 

0.066** 

0.10** 


(0.033) 

(0.034) 

(0.033) 

(0.04) 

Arson ignitions day t — 8 


0.060* 


0.012 

0.043 


(0.035) 


(0.03) 

(0.028) 

Arson ignitions day t — 9 


0.030 



0.033 


(0.025) 



(0.038) 

Arson ignitions day t — 10 


0.089*** 



0.055 


(0.03) 



(0.034) 

Arson ignitions day r — 11 





0.057 

(0.035) 

Observations 

1,915 

2,476 

2,552 

2,433 

2,413 

LL PAR(p) model 

-723.75 

-1,167.16 

-1,299.74 

-1,232.42 

-1,143.24 

LL PAR (p), all p, = 0 

-741.18 

-1,219.88 

-1,347.48 

-1,258.41 

-1,169.53 

LL null model 

-1,179.93 

-1,406.35 

-1,825.50 

-1,711.01 

-1,370.41 

Wald stat, PAR(p) 

34.87*** 

105.45*** 

95.49*** 

51.96*** 

52.58*** 

vs. all p, = 0 

Wald stat PAR(p) 

912.37*** 

478.39*** 

1,051.52*** 

957.18*** 

454.34*** 

vs. null model 







Note: Standard errors of the estimates, in parentheses, computed as the inverse of the information matrix. 

*** Indicates that Wald tests were rejected or that hypotheses tests reject the null that parameters are zero at 1% significance, ** at 5%, and * at 10%. 


magnitudes of the two sets of elasticities for 
PAR(p) and Poisson estimates are similar, so 
our estimates demonstrate, in this empirical 
case, no tendency of the Poisson model param¬ 
eter estimates to be attenuated compared to 
those produced by the PAR(p) model. 

Annual Fixed-Effects Poisson Model Results 

The fixed-effects panel Poisson count model 
estimate broadly supports hypotheses relating 
the crime of wildland arson to socioeconomic, 
ecological, and managerial factors (table 5). 
Consistent with theory and our results for 
two locations in the daily count models, 
police officers per capita are negatively and 
statistically significantly related to annual 
wildland arson counts. The poverty rate is 
positively related to wildland arson, as ex¬ 
pected. Parallel to the most recent theoretical 
developments and empirical findings for major 
crimes, wage rates are negatively related to 
wildland arson ignitions. The higher wage 
rates occurring in the late 1990s might offer 


a partial explanation for lower arson rates 
observed in our data, just as they were found 
by others to partially explain lower crime 
rates generally. Unemployment, another 
measure of the opportunity costs of crime, 
is not significantly related to arson ignitions, 
consistent with some results on major crimes 
found by Gould, Weinberg, and Mustard. The 
coefficients on the time trend and population 
are positive and negative, respectively, which 
is somewhat counterintuitive. Population's 
negative effect might be due to the loss of 
wildland areas, which usually results from 
development, where such crimes can occur—a 
loss in firesetting locations that is greater than 
the added new arsonists that an increase in 
population could bring. We cannot explain 
the positive trend in wildland arson, once all 
other included variables are accounted for. 

Consistent with previous wildfire risk mod¬ 
eling for Florida (Prestemon et al.), arson risk 
in the current year is related to fuels, fuels man¬ 
agement, and weather. Arson counts are nega¬ 
tively related to previous years’ wildfire extent, 
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Table 3. Poisson Autoregressive (Pooled) Model, Daily Wildland Arson Ignition Counts in Five 
Locations in Florida, 1994-2001, as a Function of Population-Interacted Variables 


Variables 


Variables 


Constant 

0.09 

August 

-9.38 


(0.17) 


(6.83) 

Police per capita, 

-0.027 

September 

-15.11* 


(0.034) 


(8.86) 

Poverty rate, 

2.90*** 

October 

-0.94 


(0.94) 


(6.44) 

Unemployment rate, 

3.30 

November 

-6.59 


(2.63) 


(6.61) 

Real retail wage, 

30.29*** 

Time, 

-0.55*** 


(5.67) 


(0.09) 

Wildfire years t to t — 2 

0.024 

Population, 

-0.25*** 


(0.017) 


(0.04) 

Wildfire years t — 3tot — 5 

0.027 

Arson ignitions day t — 1 

0.25*** 


(0.073) 


(0.02) 

Prescribed fire year t 

-0.31 

Arson ignitions day t — 2 

on*** 


(0.32) 


(0.02) 

Prescribed fire year t — 1 

-0.30 

Arson ignitions day t — 3 

0.07*** 


(0.34) 


(0.02) 

KBD1, 

5.43*** 

Arson ignitions day t — 4 

0.10*** 


(0.77) 


(0.02) 

Saturday 

8.06** 

Arson ignitions day t — 5 

0.026* 


(3.42) 


(0.015) 

Sunday 

2.12 

Arson ignitions day t — 6 

0.080*** 


(3.26) 


(0.017) 

Holiday 

20.84*** 

Arson ignitions day t — 7 

0.073*** 


(4.35) 


(0.017) 

January 

5.74 

Arson ignitions day t — 8 

0.027** 


(5.42) 


(0.013) 

February 

21.44*** 

Arson ignitions day I — 9 

0.042*** 


(4.87) 


(0.014) 

March 

14.58*** 

Arson ignitions day t — 10 

0.035*** 


(5.24) 


(0.014) 

April 

17.99*** 

Arson ignitions day t — 11 

0.030** 


(4.95) 


(0.014) 

May 

8.47 

Observations 

11,789 


(5.41) 

LL PAR(p) model 

-5,805.12 

June 

-9.39 

LL PAR(p), all p, = 0 

-6,125.51 


(7.64) 

LL null model 

-7,512.39 

July 

-6.08 

Wald stat, PA R (p ) vs. all p, = 0 

640.78*** 


(6.29) 

Wald stat PAR(p) vs. null model 

3,414.54*** 


Note: Standard errors of the estimates, in parentheses, computed as the inverse of the information matrix. 

***Indicates that Wald tests were rejected or that hypotheses tests reject the null that parameters are zero at 1% significance, ** at 5%, and * at 10%. 


with coefficients on six out of the first seven 
years’ lagged wildfire area negative and sig¬ 
nificantly different from zero at the 1% signifi¬ 
cance level. The twelfth year lag is positive and 
significant, indicating a returning risk many 
years after local fires. Because previous wild¬ 
fires index fuel reductions and hence higher 
costs or lower success rates of wildland arson, 
this result accords with our expectation. Haz¬ 
ard reducing prescribed fire, another measure 
of fuels, is negatively related to arson ignition 
counts, also supporting our hypotheses about 


arson ignition costs or changed ignition success 
rates. On the other hand, pulpwood harvesting 
activities, sometimes done to thin forests and 
reduce fuel loads, are positively linked to igni¬ 
tions. This can be explained by recognizing that 
harvesting can temporarily increase downed 
woody debris left over from cutting (e.g., tree 
tops, limbs), which provide fuel for fire. Arson¬ 
ists might view such areas of recent harvests 
as places where arson success could be higher, 
so they target them. Finally, the measure of 
average weather, the Nino-3 SST anomaly, is 
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Table 5. Fixed-Effects Time-Series Poisson Model Estimate of Annual Wildland Arson Ignition 
Counts in Florida, 1995-2001 


Variables 

Parameter Estimate 

Variables 

Parameter Estimate 

Police per capita T 

—0.61** 

Wildfire acres,_n 

0.20 


(0.20) 


(0.35) 

Poverty rate T 

1.75** 

Wildfire acres,_ 42 

^2*** 


(0.85) 


(0.36) 

Unemployment rate T 

-1.39 

Hazard prescribed fire acres T 

-0.84*** 


(2.86) 


(0.27) 

Real retail wage. 

-0.83*** 

Hazard prescribed fire acres,_i 

0.29 


(0.26) 


(0.26) 

Wildfire acres T _i 

-0.65*** 

Hazard prescribed fire acres T _ 2 

0.47 


(0.12) 


(0.34) 

Wildfire acres,_ 2 

-0.37*** 

Pulpwood harvests,-! 

3 20*** 


(0.10) 


(0.56) 

Wildfire acres T _ 3 

-4.72*** 

Pulpwood harvests T _ 2 

-1.10 


(0.76) 


(0.70) 

Wildfire acres,_ 4 

-2.46*** 

Pulpwood harvests T _ 3 

0.29 


(0.77) 


(0.65) 

Wildfire acres,_ 5 

_^ Q'J*** 

Nino 3 SST anomaly T 

-0.18*** 


(0-67) 


(0.04) 

Wildfire acres T _ 6 

0.21 

1998 dummy. 

0.33*** 


(0.44) 


(0.09) 

Wildfire acres T _ 7 

-1.22*** 

Time, 

0.28*** 


(0.42) 


(0.05) 

Wildfire acres T _ 8 

0.50 

Population, 

-3.05* 


(0.36) 


(1.66) 

Wildfire acres,_ 9 

-0.04 

Observations 

402 


(0.36) 

Log-likelihood, full model 

-1,129.86 

Wildfire acres T _ i0 

0.57 

Log-likelihood, null model 

-4,616.55 


(0.35) 

Wald test, full vs. null 

6,973.39*** 


***Indicates that Wald tests were rejected or that hypotheses tests reject the null that parameters are zero at 1% significance, ** at 5%, and * at 10%. 


negatively related to arson ignitions, while the 
dummy variable accounting for extreme con¬ 
ditions related to the ENSO cycle of 1997-8 
is positively related. Both of these effects were 
expected, based on previous wildfire modeling. 

Implications 

The time-series properties and general causes 
of wildland arson have not been adequately 
explored, and the research reported here, par¬ 
allel to research into crime and human-caused 
wildfires, leads to four principal findings. First, 
wildland arson demonstrates temporal cluster¬ 
ing, which supports hypotheses of either serial 
or copycat firesetting. As identified in all six 
cases analyzed, wildland arson demonstrates 
highly significant autocorrelation at the daily 
time scale, for periods lasting up to eleven 
days, just as has been shown for other types of 
crime. The apparent episodic pattern of arson 
implies that it has some short-run predictabil¬ 
ity, which should aid in developing tactical 


responses to arson outbreaks. Additionally, 
temporal clustering overlays weekly and intra¬ 
annual patterns in firesetting that could be ex¬ 
ploited by law enforcement and wildland fire 
managers to develop strategies that limit arson 
occurrence. 

Combining our first finding with our evi¬ 
dence on law enforcement leads to our second 
finding, that law enforcement resource place¬ 
ment strategies could be effective. The statis¬ 
tical results reported here add to the support 
provided by published research in the wider 
field of criminology. Agencies could preposi¬ 
tion law enforcement resources in locations 
where recent suspected wildland arson igni¬ 
tions have occurred. They could also regu¬ 
larly increase arson enforcement on days and 
months of the year when ignitions are more 
common and during droughts. Specifically in 
Florida, this means raising enforcement lev¬ 
els on Saturdays and holidays and during the 
spring fire season. Also in Florida, some de¬ 
cisions can be made months in advance by 
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monitoring forecasts of the El Nino-Southern 
Oscillation (e.g., Ji, Behringer, and Leetmaa). 

Third, we find that locations with difficult 
economic conditions have higher rates of wild¬ 
land arson, other variables held constant. The 
negative relationship with wages and positive 
relationship with poverty that we find in our 
annual model validates a contention that eco¬ 
nomic conditions matter and provides some 
evidence that arsonists behave consistent with 
an economic model of crime. Given this, an¬ 
other strategy by law enforcement would be 
to pay attention to wildland arson in times of 
economic downturns and in places with chron¬ 
ically low wages and high poverty. 

Fourth, forest management activities are re¬ 
lated to wildland arson. Sometimes, manage¬ 
ment exacerbates the risk of wildland arson, 
while other actions alleviate the risk. Although 
our statistical finding that timber harvesting 
activity is a positive risk factor is not proved 
to be causal, these results are consistent with 
other research (Prestemon et al.). In contrast, 
because fuel reductions caused by prescribed 
fire and previous wildfires are correlated with 
lower arson rates, it makes sense that managers 
could reduce arson rates by reducing available 
fuels. This finding is also consistent with an eco¬ 
nomic model of wildland arson crime, wherein 
lower fuels increase the cost of successful fire¬ 
setting. 

Our results highlight several research needs. 
First, daily models reported here are for high 
arson locations in Florida, so new research 
should explore whether the temporal cluster¬ 
ing and seasonal patterns that we identify hold 
for other parts of the country. Second, a po¬ 
tentially fruitful area of new research is on 
how to augment the received optimal wildfire 
intervention models (e.g., Gorte and Gorte; 
Rideout and Omi) to include ignition counts, 
including those from arson, and their causes. 
Wildfire costs include a fixed amount, associ¬ 
ated with ignitions, and a variable amount, as¬ 
sociated with the size, severity, and duration of 
wildfires. The inclusion of counts would there¬ 
fore capture the effects of demographic and 
socioeconomic changes over time that can af¬ 
fect aggregate wildfire risks. This is especially 
true in areas where the potential economic 
damages from wildfire are greatest (e.g., in the 
wildland-urban interface). 

Third, our models have been estimated on 
a spatial scale that could be viewed as arbi¬ 
trary, based on political boundaries rather than 
another objective metric. New research could 
evaluate the spatial scale most useful for de¬ 


tecting the effects of hypothesized influential 
factors, including economic conditions and law 
enforcement. For example, finer scale model¬ 
ing could explicitly include the locations of re¬ 
cent arson ignitions to predict future ignitions 
nearby, information more useful for targeting 
law enforcement and achieving greater arrest 
and conviction rates. 

[Received December 2003; 
accepted November 2004.] 
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